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Summary  /  Overview 

The  purpose  of  this  effort  was  to  investigate  the  behavior  of  three-dimensional 
scalar  fields  in  fully-developed,  gas-phase,  turbulent  flows;  to  further  our  under¬ 
standing  of  phenomena  that  rely  on  the  behavior  of  scalar  gradients,  such  as  aeroop- 
tical  effects,  imaging  through  a  turbulent  medium,  such  as  the  atmosphere;  as  well 
as  laser  propagation  through,  and  scattering  by,  gas-phase  turbulent  flows.  It  has 
been  a  collaborative  effort  with  the  Physics  and  Astronomy  departments  at  Cal¬ 
tech,  and  the  Jet  Propulsion  Laboratory,  and  comprised  of  several  experimental 
and  theoretical/modeling  parts  that  deal  with  scalar-fields  and  their  aerooptical 
effects  in  various  free-shear  flows.  It  has  benefited  from  equipment  and  instrumen¬ 
tation  developed  under  cosponsorship  by  AFOSR  Grant  No.  F49620-94-1-0353  and 
AFOSR/DURIP  Grant  No.  F49620-95-1-0199.  In  this  three-year  effort,  we  have 
investigated  the  structure  of  index-of-refraction  fields  in  high  Reynolds  number  gas- 
phase  turbulence,  along  with  aerooptical  phenomena  associated  with  beam  prop¬ 
agation  through  such  fields.  Flows  specifically  investigated  include  jets  and  shear 
layers.  Out  results  have  shed  light  on  a  host  of  phenomena  that  are  important  in 
the  context  of  aerooptics  that  had  not  been  previously  addressed. 
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1.  Introduction 

The  research  program,  “Whole-field  measurements  of  turbulent  flow  for  the 
study  of  aero-optical  effects,”  completed  under  AFOSR  support.  Grant  No.  F49620- 
94-1-0283,  was  focused  on  the  study  and  characterization  of  the  multidimensional 
geometrical  structure  of  scalar  fields  generated  by  fully- developed  turbulence,  with 
an  emphasis  on  investigations  of  index-of-refraction  fields  as  they  pertain  to  aeroop- 
tical  phenomena  and  the  needs  of  aerooptical  applications. 

Aerooptical  phenomena  involve  the  propagation  of  transmitted  beams,  such 
as  lasers,  imaging  beams,  as  occur  in  ground-based  astronomical  and  other  ob¬ 
servations,  observations  of  the  earth  from  space,  as  well  as  from  rapidly-moving 
air-borne  platforms.  In  the  course  of  this  propagation,  wavefront  coherence  of  the 
transmitted/received  optical  beam  is  distorted  by  index-of-refraction  fluctuations 
in  the  intervening  turbulent  medium,  often  compromising  its  utility. 

In  these  phenomena,  the  resulting  wavefront  degradation  may  loosely  be  re¬ 
garded  as  the  consequence  of  two  interaction  regimes.  These  are  the  result  of  small- 
amplitude,  relatively-long  time-scale  distortions,  accumulated  over  long  propagation 
distances,  as  occurs  in  propagation  through  the  atmosphere  (Tatarskii  1961,  Clif¬ 
ford  1978,  Roddier  1981,  Goodman  1985),  and  from  relatively-higher-amplitude, 
possibly  short-time  phase  distortions,  a  consequence  of  near-field  turbulence  near 
the  receiving/transmitting  end,  depending  on  its  characteristics  (ground-based,  air¬ 
borne,  etc.).  Implemented  and  contemplated  means  to  address  such  distortions 
depend  on  the  characteristics  of  each  type,  or  combination,  of  these  effects.  Long- 
range  atmospheric  wavefront  distortion,  resulting  in  the  scintillation  of  stars  in 
ground-based  telescopes,  for  example,  is  characterized  by  order-millisecond  time 
scales.  These  permit  the  employment  of  “rubber- mirror”  adaptive-optics  techniques 
in  a  feedback-control-loop  for  substantial  real-time  wavefront  phase-error  corrections 
(Babcock  1953,  1990,  Tyson  1991). 

Near-field  effects,  as  would  be  produced  in  traversing  through  boundary  layers 
or  shear  layers  can  be  characterized  by  time  scales  that  are  too  short  to  permit 
“rubber-mirror”  adaptive-optics  techniques.  By  way  of  example,  a  boundary  layer 
with  a  thickness  of  5bi  «  1  cm,  generated  on  the  exterior  of  an  optical  window  on  an 
airplane  moving  at  ~  250  m/s,  would  be  characterized  by  a  frequency  spectrum 
scaled  by  a  characteristic  large-scale  boundary-layer  time  of, 

tbi  =  ^  ~  4  X  10“®  s  . 

^  CO 


(la) 
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Separated  shear  layers  from  the  same  airplane  would  tend  to  be  thicker,  say,  6si  ~ 
10  —  50  cm,  with  accordingly-longer  characteristic  times,  j.e., 

tsi  =  ^  4  X  10"^  -  2  X  10"^  s  .  (lb) 

At  the  high  Reynolds  numbers  that  would  apply  in  these  scenarios,  however,  the 
frequency-content  of  index-of-refraction  fluctuation  spectra  is  very  large.  The  index- 
of-refraction  fluctuation  time  scales  extend  to  several  orders  of  magnitude  lower,  cor¬ 
responding  to  very  high  frequencies,  albeit  at  lower  amplitudes.  For  such  aeroopti- 
cal  interactions,  possible  mitigation  techniques  include  nonlinear,  phase-conjugation 
wave-mixing  optical  methods,  flow  control  to  alter  the  index-of-refraction  fields  to 
ones  characterized  by  lower  distortions,  or  to  ones  hopefully  lending  themselves  to 
easier  optical  correction/adaptation. 

In  all  these  applications,  progress  is  hampered  by  an  incomplete  knowledge 
of  the  geometrical  structure  of  the  index-of-refraction  field  in  the  turbulent  flow, 
through  which  the  optical  beam  must  propagate.  Successes  to  date,  using  “rubber- 
mirror”  optics,  for  example,  rest  more  with  the  general  power  of  closed-loop  control 
methods,  rather  than  a  good  model  for  the  system  being  controlled;  the  effects  of 
the  intervening  turbulent  medium,  in  this  case. 

The  purpose  of  the  work  undertaken  under  sponsorship  of  this  AFOSR  Grant 
and  described  below,  initially  under  Dr.  J.  McMichael  and  completed  under  Prof.  M. 
Glauser  as  program  managers,  has  been  to  investigate  the  structure  of  scalar  fields 
in  turbulent  flows  and  index-of-refraction  fields,  in  particular,  and  its  relation  to 
aerooptical  phenomena. 
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2.  Progress,  accomplishments,  and  new  findings 


2.1  Near-field  aerooptical  effects 

Weak  aerooptical  effects,  corresponding  to  optical-beam  propagation  through 
a  confined  turbulent  region,  can  be  described  in  terms  of  the  wavefront-propagation 
phase  function, 

if  =  J k(x,  t)  •  dx  ,  (2a) 


for  each  beam  (ray)  that  emerges  from  the  turbulence-containing  domain,  ft.  In 
this  expression, 

k  =  yk  ,  (2b) 


is  the  local  wavevector  of  the  beam,  with  A  the  local  value  of  the  wavelength,  as 
it  is  altered  {e.g.,  eikonal  equation)  in  direction,  k,  and  magnitude,  27r/A,  by  the 
fluctuating  index-of-refraction  field,  n(x,t). 

For  small  deflections,  we  may  neglect  the  change  in  the  propagation  direction, 
say,  along  the  x-axis,  and  approximate  this  integral,  expressed  in  terms  of  the 
transverse  coordinate,  say,  y,  at  any  one  time,  t,  i.e., 


^{y)  =  h  jn{x,y)* 


where  ko  =2t:/\q,  is  the  wavevector  magnitude  (in  vacuum). 

If  the  index  of  refraction  fluctuations  are  scaled  by  An,  e.y.,  if  they  arise  as  the 
result  of  mixing  of  two  fluids  with  an  index  of  refraction  ni  and  n2,  corresponding, 
for  example,  to  the  indices  of  refraction  of  two  freestream  fluids  entrained  in  a  tur¬ 
bulent  shear-layer  region,  then  we’ll  have  a  maximum  index-of-refraction  excursion 
of  An  =  jni  —  n2l.  If  An  is  not  too  large,  the  leading-order  beam-propagation  effect 
may  be  regarded  as  the  phase-accumulation  in  the  original  x-direction  of  propaga¬ 
tion,  say,  as  a  function  of  a  transverse  coordinate,  y.  This  will  be  scaled  by  the 
index-of-refraction  magnitude.  An,  and  the  extent  of  the  turbulent  region,  L,  in 
the  direction  of  propagation.  This  motivates  the  definition  of  a  scaled  propagation 
perturbation  phase  function. 


-/  ^  _  ¥’(y) 


t  / 
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which  can  be  used  to  characterize  weak  aerooptical  effects  in  confined,  fluctuating 
index-of-refraction  turbulent  regions.  This  will  be  used  to  illustrate  some  of  the 
aeroptical  effects  and  in  the  proposed  studies  of  near-field  interactions. 

Wavefront  functions  from  turbulent  regions  are  typically  estimated  in  the  lab¬ 
oratory  using  direct-measurement  techniques,  such  as  Shack- Hartmann  lenslet  ar¬ 
rays,  shearing  interferometers,  etc.  Quantitative  imaging  techniques,  however,  per¬ 
mit  a  direct  measurement  of  the  instantaneous  index-of-refraction  field  in  a  plane, 
i.e.,  n(x,y),  as  required  in  Eq.  3,  allowing  a  numerical  computation  of  the  scaled 
phasefront-propagation  function,  to  a  high  accuracy  and  resolution. 


Fig.  1  Rayleigh-scattering  image  sUce  in  the  mid-span  plane  of  a  low-compressibility 
(Me  =  0.15),  high  Reynolds  number  (Re  ~  2  x  10®),  shear-layer.  Gray  scale 
(right)  denotes  index  of  refraction.  Note  mixed-fluid  regions  of  near-uniform 
index  of  refraction. 
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The  results  of  sample  calculations,  performed  on  data  recorded  previously  (un¬ 
der  support  of  this  grant)  in  subsonic  (low-compressibility)  shear  layers  (convec¬ 
tive  Mach  number,  Me  ~  0.15),  as  well  as  supersonic  shear  layers  with  moderate 
(Me  Ci  0.54)  and  high  (Me  ^  0.96)  compressibility,  will  be  shown  below  (Fourguette 
et  al.  1998).  An  example  of  such  data  is  included  here  as  Fig.  1.  Exploiting  the  fact 
that  aerooptical  deflections  of  rays  are  small,  in  this  case,  these  data  were  image- 
processed  to  remove  the  streaks  caused  by  the  aerooptical  interaction,  to  estimate 
the  variable  index-of-refraction  field.  The  remaining  small-scale  streaks  in  Fig.  1 
are  a  manifestation  of  the  no-deflection  approximation.  As  noted  above,  the  index 
of  refraction  difference  in  this  flow  was  generated  through  the  use  of  nitrogen  and 
ethylene  as  the  freestream  gases. 

The  normalized  phase  integral  integration  (Eq.  3)  was  carried  across  the  imaged 
field  of  view,  with  L  pa  2S  (i.e.,  roughly,  twice  the  local  shear-layer  width)  equal  to 
the  extent  of  the  imaged  field,  and  Uoo  =  the  (low)  index  of  refraction  of  the 
high-speed  fluid.  Figure  2  depicts  the  scaled  phase  function,  ip(yfL),  for  the  data 
in  Fig.  1,  reproduced  here  to  provide  a  direct  visual  comparison. 

It  is  clear  from  the  scalar-field  image  data  and  the  resulting  phase  integral 
calculations  that  the  large-scale-structure  organization  of  the  flow  produces  locally 
positive-  and  negative-lens  effects.  In  particular,  depending  on  the  distance  be¬ 
tween  the  turbulent  region  and  the  aperture  of  the  receiving  optics,  the  effective 
positive  lenses  may  produce  caustics  along  the  way,  as  was  found  in  high-pressure 
experiments  of  turbulent  jets  using  the  same  gases  (Fourguette,  Dimotakis,  %  Ching 
1995). 

Also  shown  are  ensemble-averaged  spatial  spectra  of  the  normalized  phase- 
front-tilt  integrals, 

<“> 

i.e., 

S^iKyL)  =  J^T{<  i>iy'/L)i^[{y-y')/L]  >y.}  .  (4b) 

These  are  plotted  in  Fig.  3,  with  lines  of  increasing  solidity  indicating  increasing 
compressibility.  The  phase-front-tilt  spectra  can  be  seen  to  be  nearly  white  over 
much  of  the  resolved  spatial  wavenumber  range.  Interestingly,  compressibility  has  a 
stronger  effect  at  the  high  wavenumber  end  of  the  spectrum,  with  an  influence  that 
is  non-monotonic  with  increasing  Me.  In  particular,  the  lowest  scalar  gradients 


0.3  .  ■  ■  1  1  — ■ — I — I — » — 1 — I — 1 — I — I — I — I — > — ' — « — I 

0,0  0,2  0.4  0.6  0,8  1,0 

y/L 

Fig.  2  Subsonic  shear-layer  scalar-field  image  and  normalized  phase  integral, 

Eq.  3,  for  a  high-Re,  subsonic,  shear-layer  realization  (same  data  as  in  Fig.  1). 


(smallest  phase-front  tilts)  are  encountered  at  the  moderate  (intermediate)  com¬ 
pressibility  case.  This  is  probably  attributable  to  a  combination  of  Mach-number 
and  Reynolds-number  effects,  which  do  not  act  in  the  same  way.  Similar  behav¬ 
ior  has  been  observed  in  direct  measurements  of  molecular  mixing  in  chemically¬ 
reacting  shear  layers.* 


*  Work  on  mixing  presently  in  progress,  sponsored  by  AFOSR  Grant  No.  F49 62 0-9 8-1-0 052. 
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Fig.  3  Ensemble-averaged  spectra  of  phase-front-tilts  for  shear  layers  of  increasing 
compressibility.  Lines  of  increasing  solidity  correspond  to  convective  Mach 
numbers  of  Me  =  0.15,  0.54,  and  0.96. 

Similar  calculations  were  also  performed  on  turbulent-jet  data,  more-recently 
acquired  in  the  GALCIT  Variable-Pressure  Combustion  Facility  (VPCF),  under 
support  of  this  grant.  Pressure  in  this  facility  can  be  varied,  0.1  atm  <p<  15  atm, 
permitting  Reynolds  numbers  to  be  varied  over  a  large  range.  More  importantly, 
however,  optical  signals  can  be  enhanced  by  operating  at  elevated  pressures.  Ex¬ 
amples  from  preliminary  experiments  at  p  =  4  atm  and  p  =  10  atm  will  be  shown 
below.  These  also  rely  on  Rayleigh  scattering,  utilizing  ethylene  as  the  jet  fluid, 
discharging  in  a  quiescent  reservoir  of  nitrogen. 

The  Re  ~  3.0  x  10^  jet  data  at  p  =  10  atm  were  recorded  previously  (Fourguette 
et  al.  1995),  with  new  data  at  the  lower  Reynolds  number  as  part  of  a  more-recent 
investigation.**  Examples  are  shown  in  Fig.  4  and  record  a  slice  of  the  scalar  fleld 
and  associated  aerooptical  distortion  for  a  beam  (laser  sheet)  propagating  from  left 
to  right,  in  a  plane  containing  the  jet  axis.  The  earlier  (p  =  10  atm)  experiments 
relied  on  a  cryogenically-cooled  (1024  x  1024)-pixel  CCD  camera  (Photometries, 
Model  200  series).  The  elevated  ambient  pressures  resulted  in  higher  index-of- 
refraction  gradients  and  signal-to-noise  ratios. 


**  Jet  Reynolds  numbers  here  are  cited  using  the  ambient,  reservoir-fluid  (N2)  viscosity. 
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Fig.  4  Examples  of  Rayleigh-scattering  images  of  gas-phase  turbulent  jets  (C2H4 
injected  into  N2),  for  Re  1.5  x  10^,  at  p  w  4 atm,  1.4  <  x/d-j^t  <  7.5 
(left),  and  Re  ~  3.0  x  10^,  at  p  «  10  atm,  1.5  <  s/djet  <  6.5  (right). 
Aerooptical  effects  due  to  index-of-refraction  variations  are  evident  in  these 
images. 

The  more  recent  experiments  relied  on  a  back-illuminated,  higher-quantum- 
efficiency  (ca.  75%)  cryogenically-cooled  (low-noise)  CCD  camera  (Princeton  Instru¬ 
ments,  Model  TKB1024-1),  permitting  images  to  be  acquired  at  a  (lower)  ambient- 
tank  pressure  of  p  w  4  atm. 

As  with  the  shear-layer  experiments,  the  primary  interest  in  this  part  of  the 
proposed  effort  is  to  investigate  the  geometry  of  the  scalar  field  and  the  associated 
distortions  of  optical  wavefronts.  Note  that,  in  these  experiments,  the  fluid  Schmidt 
number  is  Sc  1  (c/.  Sc  «  2.0  X 10®  for  the  liquid-phase  jets  -  Catrakis  &  Dimotakis 
1996).  The  higher  molecular  diffusivity  in  the  gas  phase  produces  a  smoother  scalar 
field  and  lower  scalar  gradients  (at  small  scales).  In  comparison  with  liquid-phase 
data,  such  experiments  can  be  used  to  quantify  Sc  and  Re  effects  on  the  scalar  field. 

Before  an  analysis  of  scalar  level  sets  or  beam  propagation  can  be  made,  the 
aerooptical  “streaks”  evident  in  Fig.  4  must  be  removed,  as  was  necessary  for  the 
shear-layer  data  described  above  (Fig.  1).  We  have  developed  a  new  technique,  based 
on  complex  Fourier  filtering,  which  attempts  to  remove  these  streaks  by  application 
of  a  notch-type  filter  to  the  Fourier-transformed  images,  followed  by  a  subsequent 
inverse- Fourier  transformation.  An  example  of  the  results  of  this  procedure  is  shown 


Fig.  5  Examples  of  Fourier  filtering  applied  to  images  of  gas-phase  turbulent  jets 
at  Re  ~  1.5  x  10'^  and  p  ~  4  atm  (left),  and  Re  ~  3.0  x  10^  and  p  ~  10  atm 
(right),  to  remove  aero-optical  “streaks”  (c/.  Fig.  4). 


y/L 


Fig.  6  Normalized  phase-propagation  integral  computed  from  a  single,  Fourier- 
filtered  image  of  a  Re  ~  1.5  x  10^,  p  4 atm  gas-phase  turbulent  jet. 


in  Fig.  5.  While  there  are  ways  to  further  improve  the  filtering  process,  the  simple 
notch-type  filter  used  to  remove  the  streaks  in  Fig.  5  can  be  seen  to  be  effective. 
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Fig.  7  Normalized  phase- propagation  integral  computed  from  a  single,  Fourier- 
filtered  image  of  a  Re  ~  3.0  x  10^,  p  ~  10  atm  gas-phase  turbulent  jet. 

Using  the  Fourier-filtered  images,  similar  phase-propagation  integrals  can  be 
computed,  as  before  (Eq.  3).  Examples  of  such  calculations,  derived  from  the  data 
in  the  Fourier-filtered  turbulent-jet  images,  are  shown  in  Figs.  6  and  7.  For  these 
calculations,  Uoo  is  the  index  of  refraction  of  the  ambient  (reservoir)  fiuid  (N2,  in 
this  case)  and  An  =  njet  —  «oo  (c/.  Eq.  3). 

Phase-front-tilt  spectra  (Eq.  4)  can  also  be  computed  for  these  data.  They  are 
shown  in  Fig.  8  for  the  two  Fourier-filtered  images  depicted  in  Fig.  5.  Both  Re  and 
pressure  effects  can  be  expected  to  be  reflected  in  this  statistic. 

Experiments  over  a  range  of  Reynolds  numbers,  with  several  realizations  for 
each  flow  condition,  will  allow  an  investigation  of  Re  effects  on  the  scalar  fields 
and  the  resulting  beam-propagation  characteristics,  as  well  as  an  isoscalar-surface 
analysis  utilizing  the  framework  described  above,  that  was  applied  previously  to 
liquid-phase  turbulent  jets  (Catrakis  &  Dimotakis  1996). 

The  work  comprised  of  the  aerooptical  calculations  and  the  new  gas-phase 
turbulent-jet  experiments  described  above  were  performed  in  collaboration  with 
H.  Catrakis,  as  a  continuation  of  work  begun  with  D.  Fourguette,  as  part  of  the 
recently-completed  investigations,  under  this  grant.  A  more  complete  discussion  of 
these  results  is  beyond  the  purposes  of  this  proposal.  A  first  report  of  this  work  will 
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Fig.  8  Phase-front-tilt  spectra  for  ethylene/nitrogen  turbulent  jets.  Dashed  line: 
Re  ~  1.5  X  10^  and  p  ~  4  atm.  Solid  line:  Re  ~  3.0  x  10'^  and  p  c±  10  atm. 

be  presented  in  the  upcoming  AIAA  Fluid  Dynamics  meeting  (Dimotakis,  Catrakis, 
&  Fourguette  1998). 


2.2  Long-range  optical-beam  propagation  in  the  atmosphere 

Using  ground-based  astronomical  observatories  provides  a  useful  test  bed  for 
optical  beam/imaging  terrestrial  applications  from  space-borne  platforms.  While 
imaging  from  airborne  platforms  is  easier,  in  that  the  optical  sweep  caused  by  the 
high  speeds  of  orbiting  satellites  and  other  objects  do  not  have  to  be  contended 
with,  nevertheless,  both  applications  must  contend  with  the  turbulent  atmosphere. 
In  both  cases,  as  with  long-range  optical  propagation  within  the  atmosphere,  an 
important  issue  is  the  formation  of  caustics.  Caustics  are  singularities  that  represent 
a  many-to-one-to-many  transformation  in  the  propagation.  When  coupled  with 
finite  signal-to-noise  ratios,  adaptive-optics  systems  have  a  scant  chance  recovering 
beam  energy  that  has  gone  through  a  caustic. 

Examples  of  images  of  Arcturus,  used  as  a  bright  point-source  reference,  ac¬ 
quired  with  the  200"  Palomax  telescope,  with  the  Cassini-CCD  (1024^)  camera 
system  described  above,  is  depicted  in  Fig.  9.  These  were  binned  in  (2  x  2)-pixel 
regions  in  real  time  and  read  at  20franies/s  (at  a  512  x  512  pixel  resolution). 
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Fig.  9  Two  out-of-focus  images  of  Arcturus,  recorded  in  good-to-moderate  seeing 
conditions,  to  illustrate  caustic  formation,  on  the  Palomax  200"  telescope, 
50ms  apart,  cropped  to  11  x  12.8  arc-sec  (275  x  320  binned  pixels)  each. 

Figure  9  depicts  two  images  from  this  sequence  (50ms  apart),  formed  of  a 
bright  star  (Arcturus)  in  the  first  hours  of  the  night  (good-to-moderate  seeing), 
purposely  out  of  focus,  to  investigate  the  tendency  of  caustic  formation  by  the 
turbulent  atmosphere.  The  magnification  of  the  optical  system  was,  roughly,  50 
(unbinned)  pixels/arc-sec.  The  hexagonal  support  structure  of  the  200"  reflector 
can  be  discerned  imprinted  on  the  pattern  of  caustics.  Substantial  differences  in 
caustic  patterns  can  be  seen  to  have  occurred  in  the  intervening  50  ms  that  axe 
attributable  to  the  unsteady,  turbulent  atmosphere. 

This  effort,  which  began  under  this  grant,  was  suspended  because  it  soon  be¬ 
came  clear  that  imaging  limitations  imposed  by  our  Cassini-CCD-based  system 
prevented  us  from  addressing  atmospheric  scintillation  phenomena.  These  occur 
on  time  scales  in  the  range  of  10“^  s  <  t  <  10”^  s,  or  so.  We  are  hoping 

we  will  have  the  opportunity  to  resume  it,  following  the  completion  of  our  high 
framing-rate  digital-imaging  system,  presently  under  development  under  a  DURIP 
Grant  No.  F49620-95-1-0199,  “High  speed/resolution/dynamic-range  imaging  sys¬ 
tem  for  subsonic  and  supersonic  turbulent  flows,”  and  continuing  rmder  AFOSR 
Grant  No.  F49620-98-1-0052,  “Mixing,  chemical  reactions,  and  combustion  in  sub- 
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sonic  and  supersonic  turbulent  flows,”  and  NSF  Grant  No.  AST9618880,  “Ex-Post 
Facto  Diffraction-Limited  Imaging  Through  Atmospheric  Turbulence.” 


2.3  Controlled-turbulence  aerooptics  investigations 


In  an  exploratory  effort  to  develop  a  system  for  accurate  image  reconstruc¬ 
tion  from  turbulence-degraded  images,  we  have  considered  approaches  that  record 
wavefront  aberrations,  along  with  the  resulting  image.  In  these  preliminary  efforts, 
wavefront  aberrations  were  estimated  using  a  single,  out-of-focus  image.  Intensity 
variations  in  an  out-of-focus  image  correspond  to  variations  in  wavefront  curvature 
at  the  entrance  pupil.^^  In  laboratory  experiments,  we  have  explored  the  limitations 
of  making  a  single  out-of-focus  measurement,  as  well  as  a  simultaneous  in-focus 
measurement. 

Variations  in  wavefront  curvature  will  cause  light  from  different  parts  of  the 
pupil  to  come  to  a  focus  at  different  places  along  the  optical  axis.  If  i  is  the  distance 
between  the  resulting  focal  position,  as  determined  by  the  local  wavefront  curvature, 
and  the  nominal  focal  plane,  we  have, 

«  =  ,  (5) 

where  9?  is  the  wavefront  phase  function  (Eq.  2),  /  is  the  focal  length  of  the  imaging 
system,  and  k  is  the  wavenumber.  Note  that  i  =  i{x,y)  is  a  local  variable,  as 
defined  over  the  pupil  aperture-interior  coordinates.  If  the  detector  is  displaced  a 
distance  £det  from  the  nominal  focal  plane,  then, 

J_  ^  1  -  {UIIfl-‘  ^ 

where  I  =  I(x,  y)  is  the  intensity  field  in  the  out-of-focus  image  and  Jq  =  Io{x,  y)  is 
the  intensity  field  in  the  absence  of  phase  aberrations.  The  approximation  is  valid 
for  AI/Iq  =  (I  —  Io)/Io  1  and  ^det//  ■C  1.  These  equations  are  based  on  purely 
geometric  arguments,  and  are  valid  as  long  as  diffractive  effects  are  neghgible,  i.e., 
provided  ><  ^^  >  . 
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Fig.  10  Helium-jet  turbulent-flow  aberrator  setup.  Actual  images  in  ecich  case. 


To  transform  a  map  of  out-of-focus  intensity  measurements,  /(x,y),  into  wave- 
front  curvature,  we  need  to  know  /o(x,y),  the  intensity  field  in  the  absence  of 
phase  aberrations.  The  simplest  assumption  is  that  (the  mapped)  Io(x,y)  is  uni¬ 
form  across  the  pupil  or,  in  other  words,  that  there  are  no  amplitude  fluctuations. 
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Fig.  11  Simiiltaneous  in-  and  out-of-focus  imaging  system  optics. 

Note  that  static  amplitude  variations  (such  as  spatial  variations  in  the  transmission 
of  the  optical  system  or  in  detector  sensitivity)  do  not  affect  these  measurements, 
since  only  the  ratio  J/Jo  appears  in  these  equations.  In  practice,  then,  to  account 
for  static  ampUtude  variations,  Iq  can  be  tahen  to  be  constant  in  time,  with  spatial 
variations  determined  by  examining  the  intensity  with  no  aberrations  present. 

Equation  5,  relating  £  to  is  for  a  monochromatic  source.  In  our  prelimi¬ 
nary  experiments,  a  white-light  somce  was  used.  K  we  assume  that  the  phase  and 
amplitude  fluctuations  are  achromatic,  from  a  single  image  we  can  determine  (p  for 
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all  wavelengths,  by  applying  the  equations  above  at  every  value  of  k  desired. 

The  relationship  between  and  tp  was  estimated  by  nearest-neighbor  dif¬ 
ferencing,  he., 

+  ViJ+i  -  •  (7) 

Assuming  there  are  N  points  where  is  sampled,  this  yields  N  linear  equations 
with  N  variables.  Absent  additional  information,  the  transformation  from  p  to 
is  not  strictly  invertible.t  This  can  be  seen  by  noting  that  pij  =  c  (piston  term) 
for  all  {i,j)  gives  =  0  for  any  value  of  c.  Such  a  piston-term  contribution  is 
irrelevant,  however,  in  the  reconstruction. 

To  establish  a  transformation  from  '^‘^p  to  p,  a  singular-value  decomposition 
of  the  transformation  matrix  from  p  to  V'^p  was  performed,  which  allowed  a  “least- 
squares”  inversion.  Given  the  geometry  of  the  problem,  he.,  how  the  out-of-focus 
image  is  sampled,  a  single  inverse  transformation  was  used  for  every  realization  of 
the  data. 

Once  the  phase  is  known  across  the  pupil,  the  scalar  electric  field,  E  = 
with  A  assumed  constant,  can  be  treated  by  Fresnel  diffraction  instead  of  geomet¬ 
rical  optics.  To  determine  the  scalar  electric  field  at  the  out-of-focus  plane  of  ob¬ 
servation,  we  use  the  Kirchhoff  scalar  diffraction  equation  with  Neumann  boundary 
conditions  (Ref.  13,  p.  431), 

+  m)  ,  (8) 

T 

where  R  =  y/{x  -  x'^  +  {y  —  y'^  +  and  V  is  the  area  of  the  pupil.  This  cal¬ 
culation  must  be  performed  separately  for  every  value  of  k.  Knowing  the  system 
throughput  at  different  wavelengths,  the  intensity  at  each  wavelength  can  be  com¬ 
bined  to  give  the  expected  broadband  image.  The  out-of-focus  image  expected  from 
these  phase  aberrations  can  then  be  compared  to  the  out-of-focus  image  received. 
The  difference  between  the  expected  and  received  images  can,  in  principle,  be  put 
through  the  same  process  as  the  original  out-of-focus  image,  to  derive  an  improved 
estimate  of  p  at  the  pupil.  The  process  can  be  further  iterated  until  the  pupil 
phase  distribution,  when  Fresnel-diffracted  to  the  out-of-focus  plane  of  observation, 
matches  the  observed  out-of-focus  image. 


t  Boundary  information,  for  example,  would  allow  an  inversion  of  this  relation. 
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Fig.  12  Out-of  focus  (left),  reconstructed  and  actual  in-focus  images  (right),  for 
two  frames  (top/bottom). 

This  procedure  was  explored  in  a  bench-top  setup  by  imaging  a  lit  pinhole 
with  and  without  turbulent-flow  aberrations.  To  add  turbulent-flow  aberr actions, 
the  collimated  light  beam  from  the  point  light  source  was  passed  through  the  flow 


image 
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of  a  turbulent  downward- directed  helium  jet  in  air  (Fig.  10).  The  resulting  in-  and 
out-of-focus  images  were  recorded  on  our  Cassini  camera,  using  an  optical  scheme 
that  split  the  single  CCD  plane  accordingly  (Fig.  11).  Sample  data  from  two  double¬ 
image  frames  extracted  from  a  single  sequence  are  depicted  in  Fig.  12.  The  figure 
includes  the  out-of-focus  image  (left),  the  reconstructed  in-focus  image  using  the 
scheme  outlined  above  (top  right),  and  the  actual  in-focus  image  recorded  (bottom 
right).  As  can  be  seen,  the  two  are  in  qualitative  agreement. 

The  electric  field  in  the  pupil  plane  is  related  to  the  electric  field  in  the  focal 
plane  by  Fraunhofer  diffraction,  which  is  equivalent  to  a  Fourier  transform.  This 
means  that  a  Point  Spread  Function  (PSF)  can  be  determined  from  the  phase  aber¬ 
rations.  With  the  appropriate  PSF  in  hand,  an  in-focus  image  can  be  deconvolved 
to  remove  the  effects  of  the  turbulent  distortions.  In  our  case,  we  should  be  able 
to  deconvolve  the  in-focus  images  to  yield  the  point  soxuce  that  is  found  in  the 
absence  of  aberrations.  Pursuing  this  route,  we  find,  however,  that  our  PSF  does 
not  match  the  actual  in-focus  image  recorded  within  any  reasonable  errors.  The 
disagreement  between  the  estimated  PSF  and  the  measured  in-focus  image  implies 
that  a  systematic  error  has  been  made  in  the  computations. 

Some  amplitude  fluctuations  are  expected  when  observing  astronomical  ob¬ 
jects.  At  A  ~  0.7  fim,  the  atmospheric  transfer  function  due  to  amplitude  variations 
reaches  a  value  of  ~  0.9.^^  Assuming  that  all  of  the  intensity  variations  seen  are 
due  to  phase  fluctuations,  the  latter  will  be  overestimated.  We  also  note  that, 
roughly  speaking,  the  transformation  from  to  tp  in  this  geometry  amplifies 

errors  by  a  factor  of  ~  8.  Therefore,  a  modest  amplitude  fluctuations  will  produce 
relatively-large  phase  errors,  and  a  correspondingly  poor  PSF  estimation. 

The  preliminary  investigations  outlined  above  are  the  result  of  a  collaborative 
effort  with  B.  Kern  and  D.  Lang. 
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